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Abstract: We study the applicability of the Z-Sum approach to multi-loop calculations 
with massive particles in perturbative quantum field theory. We systematically analyze 
the case of one- loop scalar integrals, which represent the building blocks of any higher- loop 
calculation. We focus in particular on triangle one-loop integrals and identify strengths 
and limitations of the Z-Sum approach, extending our results to the case of one-loop box 
integrals when appropriate. We conclude with the calculation of a specific physical example: 
the calculation of heavy- flavor corrections to deep-inelastic-scattering structure functions. 
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1 Introduction 



In perturbative quantum field theory it is necessary to deal with loop diagrams when 
trying to perform higher-order corrections to scattering amplitudes. Regardless of the 
method being used, results are expressed in a basis of scalar loop integrals which need to be 
calculated. For single-loop calculations, all integrations are by now known, however such is 
not the case for multi-loop calculations. While numerical methods can sometimes be useful, 
ultimately analytical results are desirable. Furthermore, because of the growth in the 
number of loop integrals to be calculated as more complex diagrams are considered, a fully 
systematic procedure would be ideal since it would be more suitable for implementation as 
an automatic computer package. 

With that in mind, the Z-Sum approach has been proposed and is a very strong 
candidate [1, 2]. Z-Sums consist of concatenated sums with a particular structure and 
are useful because they generalize several other important functions including multiple 
poly logarithms. Furthermore, they form a Hopf algebra including several operations like 
product, convolution, and conjugation, which are instrumental in simplifying complicated 
summations. 

The method consists of two major steps. First the loop-momentum integration is 
expressed in terms of concatenated sums. Second, if the concatenated sums match a given 
pattern, they are systematically reduced to Z-Sums and ultimately multiple polylogarithms, 
in a procedure that iteratively reduces each individual sum from the innermost all the way 
to the outermost one. 

For multi-loop calculations, there are different options for performing the initial step, 
using either a Taylor or a Mellin-Barnes expansion. In the last decade, there has been a 
fair amount of interest in the method, mainly using Mellin-Barnes expansions and applying 
it successfully to specific calculations [3—24]. 

Regardless of the expansion method being used, the full multi-loop integration is per- 
formed using unphysical lower- loop integrations as building blocks, for instance three- loop 
integrations are expressed using two-loop integrations which themselves use one-loop in- 
tegrations. While in previous works a lot of emphasis has been put on the application of 
Z-Sum to individual calculations, in this work for the first time we pay particular attention 
to the procedure of systematically obtaining concatenated sums, necessary for application 
of Z-Sum algorithms, with a focus on Taylor expansions. We also perform a thorough 
survey of all physical and unphysical one-loop triangle integrations as building blocks in 
multi-loop calculations, including massive particles, with the aim of understanding how of- 
ten the method can be successfully applied, and in which way it needs to be complemented 
in order to solve a general loop integration. This allows us to identify the key missing 
algorithms in the Z-Sum machinery, paving the way for an extension of the method. We 
also for the first time present results for two-loop calculations involving massive internal 
particles using the Z-Sum reduction algorithms and Taylor series as the expansion method. 

The layout of the paper is as follow. We start in section 2 with a discussion of the 
alternative ways in which the Z-Sum algorithms can be applied to multi-loop integrals, de- 
pending on the use of either Taylor or Mellin Barnes expansions in intermediate steps. We 
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then focus on Taylor expansions and present a master formula for the method. In section 
3 we cover Z-Sum properties and algorithms, and discuss some missing steps required for 
a generalization of the method. In section 4 we develop a thorough survey of its appli- 
cability on one-loop triangle integrals used as building blocks in multi-loop calculations, 
and in section 5 we conclude by applying it to specific physically motivated calculations 
involving two-loop massive integrals necessary for the calculation of heavy-flavor correc- 
tions to deep inelastic scattering structure functions. Conclusions are presented in section 
6. Appendix A lists notations and conventions, while appendix B presents useful identities 
for manipulation of concatenated sums. 



2 Loop integrals and expansions 

In this section we discuss how to expand loop integrations in order to apply the Z-sums 
algorithms, and then focus on one particular method based on Taylor expansions. 

2.1 Alternative expansion methods 

A general multi-loop scalar integration can be written as: 

d D h d D k r 1 



(2tt) d ' ' ' (2tt) d P^ 1 ...P v s 



(2.1) 



where the Pj (i = 1, . . . , s) represent the denominators of generic loop propagators and have 
the form Pj = (K? — mf ) , with Ki a function of the integration variables kj (j = 1, . . . , r) 
and the external momenta, and mj are the masses of particles propagating in the loops. 
Moreover, we assume that both ultraviolet and infrared divergencies in the loop integrations 
are regularized in dimensional regularization, with D = 2m — 2e, where 2m is the original 
space-time dimension and e is an infinitesimal. 

In order to use the Z-Sum approach, we must express this integral in terms of con- 
catenated sums. This can be achieved by performing a denominator expansion either at 
the level of the momentum integration, where results of each momentum integration are 
expressed in terms of propagators in the following momentum integration, or at the level 
of the parameter integration, after the momentum integration has been performed. There 
are two alternative methods to perform the expansion, using either a Taylor expansion: 
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or the inverse Mellin-Barnes transformation [25, 26]: 

i i i r°° „, , (B 



T(z + c)T(-z) - . (2.3) 



(A + BY 27riA-T(c) J_ ioo v ' v ' \A 

We refer to 2.3 as an expansion, even though it is not, because it will lead to a series after 
the complex integration is performed by completing the contour over complex infinity and 
applying the residue theorem. 
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Approach A 


Approach B 


Approach C 


1 


Introduce Feynman parameters in first loop 


2 


Perform momentum integration of first loop 








Using Mellin-Barnes, 


3 


Express result of initial steps as one "artificial" 


express result as 




propagator in 


second loop 


several propagators in 
second loop 


4 


Introduce Feynman parameters in second loop 


5 


Perform momentum integration of second loop 


6 


Repeat procedure fo 


t more loops if necessary and perform all 


momentum integrations 




7 


Expand denominator 
using Taylor expansion 


Expand denominator 
using Mellin-Barnes 
expansion 


No expansion necessary 


8 


Perform all parameter integrations 


9 


No complex integration 


Perform all complex integrations 


10 


Obtain final result in terms of concatenated sums 



Table 1. Steps for multi-loop calculations using approaches A, B and C. 



As a result of the application of these expansions, the loop integral in eq. (2.1) is ex- 
pressed in terms of concatenated sums through a series of simple steps, initially expressing 
the momentum integration as a parameter integral and finally reducing the latter to beta 
functions, thus systematizing the calculation of both the momentum and parameter inte- 
grations. While the Taylor expansion may only be applied to the parameter integration, 
the Mellin-Barnes transformation can be used either at the momentum or the parameter- 
integration level, leading to three interesting distinct approaches to perform the expansion 
of multi-loop scalar integrals, which we summarize in table 1. 

Approaches A and B follow the standard procedure for multi-loop momentum inte- 
gration, and the expansion is only performed after all momentum integrations are done, 
where the factor being expanded is the denominator of the parameter integration. In or- 
der for the expansion to be legal, conditions must be met. For approach A, which uses 
Taylor expansions, (A/B) (see eq. (2.2)) will involve physical invariants and integration 
parameters, and the condition < 1 must be enforced for some choice of invariants 

(physical or not) over the whole range of integration. Conditions for approach B are only 
enforced when the complex integration is performed by completing the integration contour 
and using the residue theorem, with the result being a (concatenated) sum over residues. 
In order to complete the contour, the integrand must vanish at complex infinity, which 
will lead to the condition of convergence. 

Approaches A and B share many of the same steps and thus are very similar. All well 
defined expansions using approach A are equivalent to the cases in approach B where the 
result from the parameter integrations does not contain poles. Approach B is more general 
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however, since it also allows for cases where the result from the parameter integration 
involves poles, unlike approach A, where these cases would correspond to badly defined 
expansions. 

Approach C is different from the previous two as it introduces Mellin-Barnes transfor- 
mations before all momentum integrations have been performed, with the expansion being 
used to express the results of one momentum integration in terms of propagators of the 
next momentum integration, and not at the level of the parameter integrations. Unlike the 
previous case, Taylor expansions cannot be used in place of Mellin-Barnes transformations 
as it would always lead to badly defined expansions. 

While these methods have been successfully applied to specific calculations [3—24], a 
thorough study of each approach has not yet been performed. Such undertake would be 
necessary in order to understand the types of summations obtained from a general loop 
integration and to learn how often the Z-Sum approach is applicable. While we would like 
to eventually study all three methods in order to find out which one is most promising, in 
this work we will focus on approach A. 

2.2 General expression for approach A 

We will now derive a general expression for approach A. Following the first six steps as 
shown in table 1, Feynman parameters are introduced and all momentum integrations are 
performed. We obtain a general parameter integration of the form: 

WUdn ; (2 4) 

iu=i ^-Ef^iC( Mfc )n:=i^' tfc ' 2s - 1} (i-^)^^) d ' 

where we have performed a change of variable so that all parameter integration variables 
ui (I = 1, . . . , n) range from to 1, and we have omitted any overall factor not relevant to 
the Z-Sum procedure, even if it involves an infinitesimal part, since this is not crucial to the 
development of our discussion. In eq. (2.4), n is the number of parameter variables (equal 
to the number of propagators minus 1), 7 is the number of factors in the denominator, 
and the coefficients ]tfe ) involve invariants (dot products between external momenta and 
masses squared). The powers in the denominator terms must obey P(k,t k ,s) > 0- The 
form (1 — A) in the denominator can always be achieved by simple manipulations if not 
immediately obtained from the previous steps. Note that the form of the polynomial in 
the denominator is not unique since we are expressing it in powers of both u s and (1 — u s ) 
(s = 1, . . . ,n), with different forms leading to different (but equivalent) summations when 
expanded. 

We proceed with the Taylor expansion of the denominator followed by binomial expan- 
sions if necessary, after which all parameter integrations are of the form of beta functions. 
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All these steps can be summarized in the following master formula: 



w^n £ ^ E .... E . 



(dp 11 ^ T(i (k0) , , 

/• • • \i TT Hk,t h -i)~Hk,t k ) fj\ r (f(2«-l)) r (5(2o) 

x (t( fc> o)-t (M ), ... ,^A-D-^A.)) ! ll c (M fc ) (ll r( g(2 ,_ 1)+g(20 ) 



tfc=l 
with 

0(4 = a (e) + 1 + Yl (P(w,+l,*) ~ P(9J,^)) ' ( 2 - 6 ) 

g=i i,=o 

where P( q fl,e) = Hk,p k ) = ^ was use d to shorten the notation, and we used the multinomial 
coefficient notation: 

/ \ 1 r(ai + ... + a n + 1) 

ai, ... ,a n ! = — — — — — — — — . 2.7) 

r(ai + 1) ... T(a n + 1) 

Although eqs. (2.5) and (2.6) may not look very transparent because of their generality, 
they are extremely useful as they show the types of sums to be expected when using 
approach A. The most important part of these expressions, and what defines whether the 
Z-Sum reduction will be possible, is the argument of the gamma functions arising from the 
parameter integrations, namely eq. (2.6). It consists of a constant and a variable part. 
The latter depends on the summation variables i( q j q ) multiplied by the exponents P( q j tl e) 
of the terms in the original polynomial in the denominator of the parameter integration. 
Depending on the form of this denominator, the arguments of the gamma functions will 
involve several summation variables, and will therefore be more complicated to deal with. 
The constant term am involves integers and possibly an infinitesimal e stemming from the 
D = 2m — 2e of dimensional regularization. 

3 Gamma function expansion and Z-Sums 

Whenever an infinitesimal is present in the argument of a gamma function, we use the 
following expansions: 

n-1 

T(n + e) = 9(n > 0) T(l + e) T(n) ^ i Zj(n — 1) 

(i + O (-1 

e r(l-n) 



i=Q 
00 



+ 0(n < 0) r(1 ^ } _ ( | r £i . 



i=0 

00 



^ ' i=0 
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Conversion 



Z(n,...) <-> S(n,...) 



Limit shifting Z(n+c, . . .) O Z(n, . . .) 

Limit multiplication Z(2n, ...)■<-> Z(n, . . .) 

Multiplication -^( re > • • -)Z(n, . . .) O Z(n, . . .) 

Conjugation ^ . J ^ + ^ m S(i+rfi, . . .) -B- 5(n, . . .) 



Convolution 



- — . * mi Z{i+di, . . .)- — 2 Z(n-i+d 2 , ...)++ • • •) 

, =1 (aii+bi) (a 2 (n-i)+b 2 ) 

Convolution and conjugation 

Table 2. Properties of Z and S-Sums. 



where n is an integer, 9 is a boolean step function, and Z{ and S% are called Euler-Zagier 
and harmonic sums, respectively. These functions are special cases of Z-Sums and S-Sums 
[1, 2, 27, 28], defined by: 

Ex* 1 x ik 
-^r • • • "Jr , (3-3) 

n>ii>i 2 >...>i fc >0 1 fc 

S(n) = 0(n > 0), %mi,...,m fe ;a;i,...,^)= V . . . -Jr , (3.4) 

n>ii>i2>...>ift>l K 

where is the "depth" and {m\ + . . .+m^) is the "weight" of a sum. Z-Sums also generalize 
other functions including multiple polylogarithms of Goncharov, harmonic polylogarithms 
of Remiddi-Vermaseren, Nielsen's generalized polylogarithms, and classical polylogarithms 
[29—33]. Z-Sums satisfy a Hopf algebra, some properties of which are illustrated in table 
2. More details on these operations can be found in [1, 28]. 

Note that, at first look, the building blocks of the algorithms in table 2 match the ones 
obtained in eq. (2.5), namely binomial coefficients, x l factors, inner Z-Sums, and l +b 
factors originating from the simplification of ratios of gamma functions. The intent of the 
method is to use the properties of Z-Sums to systematically reduce the expression obtained 
from the expansion of loop integrals to multiple polylogarithms. 

3.1 General simplifications 

Before one can use the Z-Sum algorithms listed in table 2 on special cases of eq. (2.5), some 
simplifications are necessary. Some of these are trivial but are included for completeness. 
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Ratios of gamma functions After the expansion in powers of e has been performed, 
one still needs to simplify ratios of gamma functions in order to use the Z-Sum algorithms. 
This can be done using: 

T(i + a)_ 1 _ 1 fb-a-l\ 



T{i + b) (i + 6-l)(i + 6-2)...(i + a) T(b - a) ^ \ j Ji + a + j' 

(3.5) 

if (b — a) is a positive integer and: 

T( + ) a ~ b 

r (i + b) ={i + a - W + a-2)...(i + b) =J2^a-b-j(.b, ... ,a - 1) , (3.6) 

if (6 — a) is a negative integer, where we define the kappa function K s (t) as the sum of 
products of every subset of t with s elements, for example: 

Kq(. • •) = 1 , 

Ki(a, 6, c) = a + b + c , 

K2(a,b, c) = (a& + ac + 6c) , (3-7) 
K3(a, 6, c) = abc . 

Partial Fractioning If more than one pair of gamma functions exist, one needs to 
perform partial fraction using: 



(ai + b) mi (ci + d) m * 

ad-cb) m 2 -l ) \ad - cb ) (ai + b)i K '' 

— 'mi + 771,2 — J — 1 \ / a \ 1 



+ ( nr/ - rh I ^ 



ad — cb J ^ \ mi — 1 I \ad — cb I (ci + d)i ' 
i=i 

where mi, rri2, a, and c are integers. 

Lowering and Raising Operators In some cases we need to perform sums involving 
powers of the summation variable in the numerator. Sometimes it is useful to remove these 
factors by using derivatives, as in: 

E^/« = (^J (5> i/(i) J • (3 - 9) 

Similarly, when the i m factor is in the denominator we could use: 
i=l " ' \i=l / 
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In case f(i) also depends on x we can always make the substitution x % — > (xqx) 1 (but not 
in f(i)), perform the derivative or integration with respect to xq and ultimately take the 
limit xq — > 1. These are called "lowering" and "raising" operators [2]: 

(aT) • f(x) = x ^f( x ) , 
f' x dr' 

(*+) • f{x) = / ^/(*') , (3.11) 
Jo x 

(x+r-l = -^ln m (x) . 

Reversing summation order If a sum involves the summation variable i solely as 
(n — i), where n is the upper limit, it can be trivially simplified as: 

n n 



£/(n-t) =£/(<) • (3.12) 



i=0 i=0 



Shifted sums The application of step functions 8 on sums will modify the summation 
limits. In order to use the Z-Sum algorithms, these need to be shifted back to the standard 
values. Some examples of such an operation are: 



£ 9(i < a)f(i) = 6(n < a) ]T /(*) + 0(0 < a < n) /(») , 

i=0 i=0 i=0 

n n a 

Y,0(i>n-a)f{i) = e(n<a)Y J fii)+0(n>a)Y,f{i + n-a) , (3.13) 

i=0 i=0 i=0 

n n n—a 

£e(*>o)/(i) = e(a<0)Y,f^)+G(0<a<n)Y,f(i + a) . 

i=0 i=0 i=0 

Binomial synchronization After sum shifting, the binomial coefficient, if present, might 
not match the summation's limit. It can be "synchronized" using: 

vV n W-uM-vV n + a> \ r(n + 1) r(i + 1) r(n-i + a + l) 

frfV + bJ A J_2 ^V i T(n + a + l)T(i + b+l)r(n-i-b+l) ni + a) ' 

(3.14) 
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3.2 Reduction to Z-Sums algorithms 

Combining the Z-Sum algorithms of table 2 with the simplifications of subsection 3.1, we 
obtain reduction algorithms for the following summation types: 



n 



^ (aii+6i) m i (a r i+b r 
i=i 



X(c 1 i+d 1 ;...)...X(c s i+d s ;...) , (3.15) 



n ■ f n 

B E ^ff^*™ — - ^ X{c 1 f?+d 1 ;...)...X{c s f?+d s ;...) , (3.16) 

i 



n 



C V h A— ■ ■ ■ / i\ X( Cl i+di; ...)••• X{c s i+d s ; . . .) , (3.17) 

where each X represents either a Z-Sum or S-Sum, 6^ and di are bounded numbers, 1 and 
/" equals either i or (n—i), with at least one of each present for algorithms B and D. These 
algorithms are similar to those presented in [2] with the addition of ai and q multipliers, 2 
and are used iteratively for a systematic reduction of sums obtained from loop integrals, 
applying them from the innermost sum all the way to the outermost one. We will refer to 
expressions where algorithms A, B, C, and D are sufficient as having a seamless reduction. 

Missing Algorithms While algorithms A, B, C, and D allow for the solution of many 
loop integrals, a study of eq. (2.5) shows that they are not sufficient for the general case 
when using Taylor expansion. The problem stems from two main limitations, namely a lack 
of a multiplier a>£ in algorithms C and D for some parameter integrations with quadratic 
denominator, and the requirement that the offsets bi and di be bounded numbers. 

If we focus on one limitation at a time, it is possible to list the missing algorithms as 
given by: 

Higher Power ^ n M x Z(2i,...) , (3.19) 

Variations of Type C i=o V 1 J [ai + > 



ni 



Unbounded Shifting of ^ - — -Z(i - 1, . . .)Z(i + n 2 , . . .) , (3.20) 
Z-Sum Upper Limit i=i 1 ^ 



ni 



Unbounded Offset 
of Type A 



V - . Z(t-1;...) , (3.21) 
^ (i + n 2 ) m 



i 
m 



Unbounded Offset 

of Type C i=i 



y 71 * . S(i;...) , (3.22) 



1 For the purpose of this paper, we will use the term bounded numbers to define numbers that involve 
integers and variables from summations with a finite number of terms, while we will denote as unbounded 
numbers that involve summation variables going to infinity. 

2 The algorithms in eqs. (3.15)-(3.18), without a; and c; multipliers, have been implemented in a very 
useful FORM [34, 35] package presented in [36]. 
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where b is a bounded number while n\ and 712 are unbounded, and we highlighted in red 
bold face the stumbling block in each case. 

In eq. (3.19), we have used a general factor a in the denominator (ai + b) m . This factor 
comes straight from the powers of u s and (1 — u s ) in the denominator of the parameter 
integration (see eq. (2.4)). For the cases involving single-loop triangle integrations, to be 
discussed in section 4, only solutions with a = 2 are necessary, since those polynomials are 
at most quadratic. 

It is important to note that more missing algorithms exist when combining more than 
one stumbling block. While some individual calculations involving sums like the ones 
above can be expressed in terms of multiple polylogarithms when all concatenated sums 
are taken into account, 3 a systematic and iterative procedure to reduce these types of sums 
to Z-Sums would be very useful but has not yet been developed, and in fact it is unknown 
whether it is possible. 

4 Triangle loop integrals 

In order to study the applicability of approach A we will consider the case of all one-loop 
triangle integrals. They are the simplest yet non-trivial class of loop integrals and are 
important since they are building blocks of multi-loop calculations. In order to obtain the 
results presented in this section we have developed extensive computer algorithms using 
Mathematica as a programming language. While the results discussed refer specifically to 
triangle integrations, the same procedure is applicable to more complicated integrations 
like box integrations and beyond. 

4.1 Reducible loop-like integrals 

We have used two different methods for evaluating triangle loop integrals, which we call 
"forward" and "backward" approaches. 

Forward approach When faced with a specific triangle loop integration, one possible 
way to find out whether the Z-Sum method would be sufficient is to use what we call 
forward method. In this approach, we consider every possible parametrization of the 
integration, study different ways of performing the expansion, and verify whether one of 
the resulting sums matches the Z-Sum algorithms, thus leading to an analytic expression 
in terms of multiple polylogarithms. In general, each diagram will lead to one or more well 
defined expansions, some of which may immediately match the structure of the Z-Sums 
algorithms, thus leading to a seamless reduction. For other cases, however, none of the 
expansions obtained can be immediately reduced. If the result is known by some other 
integration method, it is always possible to manipulate the expression at the summation 
level, using the algorithms presented in Appendix B, to obtain the known result. The 
steps for performing such manipulation, however, are always different for each case, and 
a systematic generalization of the procedure, necessary for obtaining new solutions, is still 
lacking. 

3 Appendix B lists several operations useful when performing operations at the summation level. 
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Sum 


Beta Functions 


V^oo r(i+d) i 

Z^i=o r(i+i)r(d) x 


[0,0], [i,0], [2i,0] 


EIU (?) 


[0,0], [i,0], [n-i,0], [i,n-i] 


E"=o /(M) x * 


[0,0], [i,0], [2i,0], [n-i,0], [2n-2i,0] 


EU ( q 7) 


[0,0], [i,0], [g-n-i,0], [i,g-n-i] 


E£ n f(q- n >i) xi 


[0,0], [i,0], [2»,0], [g-n-* s 0], [2g - 2n - 2*, 0] 


e£ ( 2 d 


[0,0], [«,0], [2n-»,0], [i,2n-i| 


ESo /(2".0 


[0,0], [i,0], [2«,0], [2ra-i,0], [4t»-2*,0] 




[0,0], [i,0], [2g-2n-i,0], [i,2g-2n-t] 


E'lo 2 " /(2g-2n,i) x< 


[0,0], [i,0], [2i,0], [2g-2n-i,0], [4g-4n-2*,0] 



Table 3. Single sums used as building blocks, obtained from Z-Sum algorithms A, B, C, and D 
(eqs. (3.15-3.18)). Each bracket [pi,p2] represents a beta function (see eq. (4.2)), and so the correct 
number of bracket needs to be selected corresponding to the number of parameter integrations (2 
for triangles, 3 for boxes, etc). See eqs. (4.1) and (4.3) for the relation between x and [pi,P2]- In 
the last two rows q represents the upper limit in the hidden n sum. f(j,i) stands for any of the 

following: 1, rffffw ' r(H+'i)t) or ru-i+i)r\d) r J+ffw ' where d wiU become the denominator 
power in the loop-like integration. 



The forward method is interesting because it gives the form of summations that can- 
not be seamlessly reduced, and has been key in understanding what types of necessary 
algorithms are still missing (eqs. (3.19-3.22)), however it is not very convenient for find- 
ing all possible summations that can be seamlessly reduced since it requires a brute force 
approach. 

Backward approach Another more systematic approach for finding all possible solu- 
tions is to consider the Z-Sum algorithms A, B, C, and D as building blocks, and use them 
to construct all possible concatenated sums that can be seamlessly reduced. We then per- 
form the calculation in reverse order, obtaining in the end all loop-like integrations with 
systematic reduction. We start this procedure by listing all building blocks, as shown in 
table 3. 

Each variable x l in table 3 has one of the following forms: 

C\ C l \p 1 ,p 2 ], C l [p 1 ,p 2 ][p3,P4], C t [pi,p 2 ][p3,P4}--- , (4.1) 

where c* is a coefficient, while \pj,pk] can be chosen from the corresponding line in the 
second column of table 3 and represent a beta function defined by: 

= T(p + a )T(p k + a k ) = ^ (1 _ n)Pfe+afc - 1 ^ 

T(pj +p k + a jk ) Jo 
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with aj and constants. If we consider the product of all x variables chosen, we obtain: 

xl 1 ... X \f = 4 1 . . . C\f X [pi , P2 ] . . . [p26_i, p 2 b] , (4.3) 

where the subscripts in x l and c l differentiate between the individual sums picked as 
building blocks, s is the total number of sums, b is the desired number of parameter 
integrations, and the product [pi,p 2 ] • • • [P2b-i,P2b] consists of a combination of any of the 
brackets allowed in each building block. 

While there are other building-block sums that lead to a seamless reduction, for exam- 
ple: 

12 1 I « I I I I 2=» I 

L 2 J /\n\\ L 2 J L 2 J / \ q— n\\ L 2 J 

E( L ?V" E/diJ.-V. E ( L TK E m^J.-V, Mi 

i=0 ^ ' i=0 i=0 ^ ' i=0 

they do not lead to a desirable integrand and so are not included in table 3. 

In order to make the approach more transparent, let us consider the following example. 
Pick two individual sums: 

V , r(<+ jL *i and T (') 4 , (4.5) 

^ r(i + i)r(d) 1 Vj7 2 1 J 

and choose a combination of allowed [pi,p 2 ] representing beta functions, for example [2i, 0] 
and [j, 0] (two since we are discussing triangle integrations), as prescribed in table 3. Con- 
catenating the two chosen building blocks we obtain: 

T(i + d) ^(21 + 01 + 1)^02 + 1) y> /A j r(i + fl3 + i)r(a 4 + i) 

^ r(i + i)r(d) Cl r(2i + o 12 + 2) ^ [jj 2 r(j + o 34 + 2) 



1 du uf(l -ni) a2 U2 3 (l - u 2 ) a4 
(1 - ci uf (1 + c 2 n 2 )) a! 



(4.6) 



which is a loop-like integration and a special case of eq. (2.1). Note that since there are 
enough degrees of freedom (c± and c 2 ) the expansion is convergent. 

Going back to the building blocks presented in table 3, we now would like to consider 
every possible way to concatenate them. A very useful method is to visualize them using 
rooted trees as a guideline. 4 A rooted tree is composed of knots and links between knots. 
A connected rooted tree has a single root on top. A knot may be linked to several lower 
knots, but only one upper one. 



4 A related use of rooted trees and their relation to Hopf algebras can be found in Appendix A of [1]. 
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U 


I 


(t) (2i) 


A 


(i,2i) (2*, 2i) 






(i,2j) (2i,j) (2i,2j) 


A\ 


(i,i,2i) (i,2i,2i) (2i,2i,2i) 


/ 


\ 


(i,i,k) (i,i,2k) (i,2i,k) (i,2i,2k) [2i,i,k) (2i,i,2k) 
(2i,2i,k) (2i,2i,2k) 


; 


\ 


(i,j,2i-2j) (2i,j,j) (2i,j,2i-j) 

(2i,j,2j) (2i,j,4i-2j) 






(i,j,k) (i,j,2k) (i,2j,k) (i,2j,2k) (2i,j,k) (2i,j,2k) 
(2i,2j,k) (2i,2j,2k) 



Table 4. All fully concatenated summations (up to four sums) using connected rooted trees as a 
guideline. The terms in parentheses represent the upper limit on summations from second to last, 
with the outer summation going to infinity. For example: (i, 2j, k) = J^iLo Sj'=o SfcLo Ez=o- 



#Sums 


Rooted Trees 


1 




2 


I 


3 


' I A 


4 


I 'A *| II A\ A I 



Table 5. All rooted trees up to four knots built by joining connected rooted trees. 

We will be considering rooted trees involving up to four knots, since this will be 
sufficient for the present discussion, however we note that rooted trees can be related to 
summations with any number of concatenated sums. We start by finding all connected 
rooted trees including up to four knots, and from these obtain all concatenated sums, as 
shown in table 4. 

Next we need to also consider trees with disconnected branches, equivalent to summa- 
tions involving more than one sum with infinity as the upper limit. This can be done by 
combining connected trees with the correct number of knots, as shown in table 5. 

Following this procedure, we obtain the results presented in table 6. These results 
may be used to compare with the parameter integration obtained from specific diagrams. 
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Table 6: Parameter integrations leading to summations with seamless reduction. 

Ine integration torm is given by: J ! 1 jp 1 , wnere 

D is the denominator presented in each row and ae are constants. represents 

££=o#, where c e is a given coefficient, dt = r(d)rffi) ' % = ^^^grj , 
and 6j = 6°, with s equal to 1 or 3. 



# Denominator (D) 



Sum 



1(1- Citt 2 ) 



EiT^&o&i 



2 (1 



ciu|) 



J2u d h b ob2h 



3(1- C1U1M2) 



C1U1U2) 



Ei? ' dnhMh 



5 (l - c x u\u\) 



6 ((1 - ciiti) (1 - c 2 u 2 )) 



7 ({I - c lUl ) (1 - c 2 up) 



((l-ci«?) (l-caul)) 



E? Yli^ di 1 di 2 b 2 i 1 b2i2 



9 ((1 - ciit 2 ) (1 - C1C2M2)) 



E? Ej 2 di 1 -i 2 di 2 bob il 



10 ((1 - C1U1U2) (1 - C1C2U1M2)) 



Ei? dj 1 -i 2 di 2 bi 1 bj 1 



11 ((1 - Cilti) (1 - C1C2W1U2)) 



Eii di 1 ~i 2 di 2 bi 1 bi 2 



12 ((1 - ciui) (l - cic 2 uitti)) 



ii Z— ri 



dii —i 2 d{ 2 Ui 1 02i 2 



13 1 


Jl-cil 


Ml £2Ml 1 

V "2 y 




1 Ei^Eia ©dii&ii-iiAa 


14 1 


1 - Ci 1 


/ -1 , C 2 Ml (1—112) ^ 


i u 2 j 


1 Eii Ei 2 {i 2 )diibi 1 _i 2 bi 2 



15 (1 - ci (1 + C2U1) W2) 



E~EjLfiKM 



12 



16 (l - Ci (1 + C 2 Ml) u| 



E?EI 2 (ij Hi^iA: 



17 




ClUl (l+ C2(1 ^^)) 


Eii Ei2 {il)diibi 1 -i 2 ,i 2 b2i 2 


18 


p 


d«i(i + a S^ i )^) 


Ei?Ei 2 (i 2 Hi fo 2ii-2i 2 ^i-i 2 


19 


( 1 - 


ci (i + ^) 2 ui) 


Ei^Eia 1 Ci^diMh-hb^ 



20 1 - ci 1 



Eh Eg 1 



*2 A . 

2ii-i 2 w *2 



21 ( 1 — Cl (1 + C 2 Ml) U 2 



E?E? 



22 (i-ci(i+c2«ir«i 



23 



E?E? C£K«*i* 



12 



(1 - jjjjj (1 + ^f 1 ) 2 j) Ei? Eg 1 (f Rfr 

24 ((1 - ciiti) (1 - C1C3U1) (1 - C2M2)) 



5 2ii -i2°2ii -«2 



Z-^ii Z^i2 Z-<i 3 "-12 "11-13 "^3 ^2 



25 ((1 - crui) (1 - C1C3M1) (1 - c 2 m^)) 



E? E? Eii d^di^di^b 



h u 2i 2 



26 (l-Ci(l + C2Ui)(l + C3«2)) 



A3 



27 (l 


- ci (1 + c 2 Ui) (1 + c 3 w 2 ) 2 ) 


V-VOO 

Z_/ii 


Z^i 2 


E? ©(tHAAa 


28 (l 


- Cl (1 + C 2 Ml) 2 (1 + C 3 W 2 ) 2 ) 


V-VOO 

2^i! 


v^2i 
Z—'i 2 




29 ((1 


- Cllt 2 ) (1 - ClC 2 (1 + C3U1) U 2 )) 


V-VOO 

Z_/ii 


y^il 
Z^i 2 


Ei3 (13)^11 -12 ^2^1^*3 


Continued on next page 



Table 6 continued 



# Denominator (D) Sum 



30 | 


[(l-c lC2 u 2 ) (l-ci(l+^ 


) j 


\ y^oo y^ii (ii-l2\j. . J. h. h. ■ ■ 


31 | 


(1 - CiU 2 ) 


fl - cic 2 (1 + c 3 u 1 ) 2 u 2 ) 




32 | 


v l-d(H- 


- c 2 ui (1 + c 3 u 2 )) 2 


) Eij E^ 1 Ei3 ( ia 1 ) (43)^*1^2^3 


33 ! 


fl-ciu?( 


i I c 2 (l+c 3 u 2 ) \ 2 | 


Si! E^ 1 (it) Qg)^il^2il-»2^3 


34 


U-cml ( 




h) 


v^oo v^2ii v^*2 /SiA /i2\ j 1/2 /. 

/Lyj 2 Z^i 3 \ i 2 ) \i 3 ) u li u 2i 1 -i 2 Ul 2-ii 


35 (l-ciui) (1 


-C1C3U1) (1-C 2 -U 2 


) (1-C2C4U2) Er°T E 4 ^ E£ E2 d^-i^d^-i^bi^ 


36 | 


(l -<*(! + 


■ca(l + C3«i)(l + 


C4M2)) 2 


) Eij Eia 1 Ei3 Ei4 (ia 1 ) G3) O^ii h 3 h 4 



37 (I-C1C2U2) (l-ci (I+C4W1) (l + i+ff^-) 

v^oo v^ii >r^ji— 12 V^«i-i2— «3 /il— «2\ fil— i 2 — 13\ J A h h 
Z-Oi Z^i 2 £-^3 ^4 I jo )\ U / 1 *2 *2 z l * 



The denominators presented in table 6 are the result of applying the backward ap- 
proach with two beta functions, and have been obtained with a software package we 
developed using Mathematica. They represent all loop-like integrations involving three 
propagators that may be reduced seamlessly using the available Z-Sum algorithms A, B, 
C, and D (eqs. (3.15-3.18)), although not all will necessarily correspond to a loop diagram. 
No new solutions are found involving more than four sums. 

4.2 Survey of triangles in QFT 

Now that we known which loop-like integrals can be seamlessly reduced, we would like to 
cross check this result with actual loop calculations. In order to consider every possible 
triangle in QFT we assign all possible values to the six invariants (p\,p\,p\,m\,m\,m^), 
where pe represents the external momenta and m are the masses of internal particles. If we 
consider only loops with at least one external leg off-shell, this gives a total of 101 triangles, 
out of which 72 are present in the Standard Model, as shown in table 7. 

The results presented in table 7 show that only a fraction of the total number of trian- 
gles (18) can be reduced to Z-Sums using Taylor expansions and the reduction algorithms 
currently known (eqs. (3.15-3.18)). All other diagrams lead to well defined expansions that 
cannot be systematically reduced to Z-Sums without solutions being found for the missing 
algorithms (eqs. (3.19-3.22)). 

For some calculations where the Z-Sum approach with Taylor expansion could not be 
successfully applied, results may be known by other integration methods. These results can 
be related at the summation level, however the steps required for making the connection are 
not systematic, changing for every diagram and every order in e, and so are not considered 
seamless reductions because they defeat the purpose of developing a systematic approach. 
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Table 7. One-loop triangle diagrams in QFT ordered by degrees of freedom, where duplicates 
due to rotation or reflection symmetries have been excluded. Bracket refer to the invariants 
{p\,p\,p\,m\,m\,m\), the external momenta and internal masses. Square brackets (red) repre- 
sent cases present in the Standard Model. The third subcolumn (labeled S) shows diagrams with 
seamless reduction, with the number indicating the corresponding polynomial in table 6, and 
indicating no expansion necessary. 

4.3 Boxes and beyond 

The procedure presented in this section can be directly applied to parameter integrations 
originating from momentum integrations with more propagators, like boxes or multi-loop 
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integrations. The difference is that they would involve a larger number of parameter- 
integration variables, and so the denominator would potentially involve a larger number of 
terms. Both the forward and the backward approaches are applicable, where in the latter 
we would use a larger number of beta functions in order to find the results equivalent to 
table 6. While we could proceed with this calculation, we do not believe it would be very 
insightful, since we expect similar results to those obtained for triangle loops. Specifically, 
we believe some of the calculations would be feasible using the existing algorithms A, B, C, 
and D (eqs. (3.15-3.18)), while other would require a solution for the missing algorithms 
discussed previously (eqs. (3.19-3.22)). 

5 Application 

In this section we would like to focus on the successes of the approach by applying it 
to a problem of physical interest, and for the first time present results for multi-loop 
calculations involving massive internal particles using the Z-Sum reduction algorithms 
and Taylor expansions. One interesting application is the calculation of the heavy-flavor 
corrections to Deep Inelastic Scattering (DIS) structure functions, which can be relevant 
(up to order 20-40%) for small values of the Bjorken variable x [7—14, 37]. 

The scattering amplitude can be represented on a basis of matrix elements of universal 
operators whose coefficients are process dependent and can be calculated perturbatively in 
QCD. The effects of massive quarks factorize entirely in the operator matrix elements [37] 
and has been calculated analytically in the asymptotic region Q 2 ^> m 2 at 0(a 2 ) [8—13], 
while partial contributions have been obtained at 0(af ) [14]. In both cases, the effects are 
sizable and definitely to be included in the theoretical predictions to be safely compared 
with the current experimental results. These calculations include the most important set 
of contributions from operators with external heavy-flavor and gluon lines and have been 
performed using inverse Mellin-Barnes transformations (approach C, table 1). As a test 
of the applicability and potential of the Taylor expansion method we have reproduced the 
results presented in [8] . 

In practice we will have to calculate two-loop integrals with massive internal particles 
and operator insertions with Feynman rules given in figure 1 (in order to compare with 
[8] we adopt Mellin-space notation). We will start the calculation with two pedagogical 
examples that do not include operator insertions, as shown in figure 2. These diagrams are 
similar to self-energy corrections to a massless particle, however the external momentum 
obeys the identity p 2 = 0. We then proceed to the calculation of diagrams with operator 
insertions, involved in the 0(a 2 ) corrections to the operator matrix elements. These are 
separated in two groups, those that do not involve a bubble loop as a building block 
and those that do, as shown in figures 3 and 4. We will focus our attention on the first 
group since these lead to the most complicated integrals. Regarding the second group, the 
presence of a bubble loop generally leads to simpler calculations and so we will not be 
performing them explicitly, although we did verify that the method may be applied in all 
cases. These diagrams represent only the integration part of the calculation, and so we do 
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Figure 1. Feynman rules for operator insertions, where A is a light-like vector and N is the 
Mellin-space variable. 




not specify which particles appear in each diagram, but only whether they are massive or 
not. 



5.1 Diagrams without operator insertions 

In this subsection we will calculate the integrations originating from diagrams shown in 
figure 2. These calculations do not include operator insertions and are similar to self- 
energy corrections to a massless particle, however with external momentum p 2 = 0. This 
subsection is included as a stepping stone in understanding the integrations with operator 
insertions to be performed in subsection 5.2. 

For clarity, we list here the steps we will be following: 
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A-I B-II B-III 




Figure 3. Comprehensive list of two-triangle loop diagrams with operator insertion stemming from 
the calculation of the heavy-quark effects in the renormalized operator matrix elements of DIS at 
0{al). 

• perform first momentum integration and express result in terms of an artificial 
propagator in second integration; 

• perform second momentum integration; 

• expand parameter integrand using Taylor series and perform integration; 

• expand gamma functions in powers of e and apply Z-Sum reduction algorithms; 

• obtain result in terms of polylogarithms. 

5.1.1 Diagram A 

We start with an explicit calculation of the integration for diagram A (figure 2). Since we 
are interested only in the integration, we ignore all prefactors, which potentially include 
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Figure 4. Remaining two-loop diagrams necessary for the calculation of the heavy-quark effects 
in the renormalized matrix elements of DIS at O(a^). Because of the presence of a bubble loop, 
these calculations are generally simpler than those presented in figure 3. 



coupling constants, color factors, and invariants. The initial expression is given by: 
d D k 2 1 



(2vr) D (q-m^Y 4 {{k 2 +p) 2 -m? 
d D h 1 

(2vr) D [kl-w?) Vl ([k l + k 2 fy 2 ((A* - p) 2 - m 2 ^" 



(5.1) 
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perform the k\ integration and re-express the result in terms of an artificial propagator: 
i( _l)^23 r (ma- f) f 1 dx 5 123 x?- 1 x?- 1 x^- 1 

(4tt)^ r(i/i) r(i/ 2 ) r(i/ 3 ) J Q {- X2 (i - x 2 )) ui23 ~^ 



x l - w , (5.2 



(2tt 



(A; 2 - m 2 )' 4 ((Ms, +p) 2 - m 2 ) 5 (fe 2 + ££V - - 



where <5i 2 3 = (5(1 — x\ — x 2 — X3) and dx represents the integration on all xg (£ = 1, 2, 3). 
At this point the k 2 integration may be performed to obtain: 

(-ir i2345+1 i> 123 45 - p) f5 31 

(4ir) D (m 2 ) Ul2345 - D T(vi) r(i^) I> 3 ) r(z/ 4 ) I> 5 ) 

11 6 123 #456 x'j' 1-1 Xg 240 2 1 (1-X2) _i/123 + T ^3 3_1 a^ 4 " 1 Z^ 5 " 1 X^ 23 2 

(x 2 (l-x 6 ) + x 6 ) 1/12345 - D 



The integration delta functions are removed by applying a change of parameter variables: 



-ir i2345+1 r(v 123 45-D) f 1 J -, 2+ f-l 



V 2 1 (1-ui)^ 245 -^- 1 







(4vr) D (m 2 ) 1 ™^ r(i/i) r(iAj) I> 3 ) T{v A ) 1> 5 ) 
x n^- 1 (1-^-1 U3- 4 - 1 (l-ng)^- 1 V"- 1 (l-n 4 )^-f-i ( 

(l-«i n 4 )^ 12345 ^ 

The denominator may be expanded after the inclusion of a regulator a, as in (1 — u\U^) — > 
(1 — auiUi), leading to: 

( _ ir2345 +i r> 245 - f ) i> 123 - f ) 



(4vr) D (m 2 )" 12345 ~ D r(i*j) r(i^ 3 ) J> 45 ) 
^ i r(i + ^12345 - £>) r(i-^ + f) 



i=0 



r(* + i) r(* + 1^2345 



Equation (5.5) may be expressed in terms of Z-Sums and ultimately polylogarithms even 
for generic values of (£ = 1, . . . , 5) and D by using identities 3.5 and 3.6. However for 
simplicity we will finish this calculation using the specific values of U£ = 1 (I = 1, . . . , 5) 
and D = 4 — 2e. The expression in eq.(5.5) becomes: 



(4tt 



) 4 m 2 \m 2 } • ' r(i + l) T(i + 3 + E) 



It is interesting to notice that this expression is simpler than the one obtained when using 
Mellin-Barnes splitting [8], as it involves a single infinite summation, with both leading to 
the same result. 
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Since eq. (5.6) is finite in e, we apply the limit e — > 0, leading to: 

i « r(i + i) 



(47T 


) 4 


to 2 








(47T 


) 4 


to 2 




t— ' 




(47T 


) 4 


to 2 



i=0 
oo 



r(» + 3) 



Va l ( — — ) (5.7) 



i=0 



1 / x M 1 
- - log (1 - a) I 2 

The final step is to remove the regulator a. In the limit a — > 1 the last terms vanish and 
we obtain: 

T A k 2 ■ ( 5 - 8 ) 
[Air) m z 

The simplicity of the result is due to the fact that this is a two point function with massless 
on-shell external particle. 

5.1.2 Diagram B 

The calculation for diagram B (figure 2) follows the same procedure and so we omit the 
intermediate equations. The initial expression is given by: 

d D k 2 1 



(2tt) d (fc 2 - m 2)-4 ^ k2+p) 
d D k\ 1 



2 2 



(5.9) 



"3 



(2tt) D ^ h + k2 f _ m2 y 2 ^ h _ p y 

We proceed by performing the momentum integrations, expanding the denominator after 
inclusion of a regulator a, and finally performing the parameter integrations. We obtain: 

(-1)^345+1 r(*i34s - # ) r(^i 23 - §) 



(4tt) D (m 2 )" 12345 " D T(u 2 ) I> 13 ) i> 45 ) 

^ r(» + ^12345 — -D) r(i - U 13 + f ) 



(5.10) 



r(» + i) r(i + ^2345 - f ) 

After setting vi = \ {I = 1, . . . , 5) and I? = 4 — 2e the expression in eq. (5.10) becomes: 

r(2 + e) r(i + e) ^ . r(i + 1 + 2g) r(j - g) 
(47r) 4 - 2£ (m 2 ) 1+2e ^ r(* + i) r(i + 3 + e) ' l °- j 

Unlike diagram A, this expression is not finite so we need to expand the gamma functions 
in powers of e. This can be done using the method discussed in section 3, after which we 
obtain: 

r(2 + £ )r(i + g ) / z i 3_j__j_/ 2 i 

(4vr) 4 " 2£ (m 2 ) 1+2£ V 2e 2 2a 2a \ a a 1 



+ o----(!-r + -2 )Iog(l-a)) • (5-12) 



When we take the limit a — >• 1 the last term vanishes and the final result becomes: 

^ + e)r(l + e) ( 1 + l) . (5.13) 



(4vr) 4 - 2£ (m 2 ) 1+2e V 2e 
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5.2 Diagrams with operator insertions 

In this subsection we will reproduce results presented in [8] by calculating diagrams in- 
volving operator insertions, shown in figure 3, first obtained by [37]. We perform the 
calculation explicitly for two examples, and also present results for the remaining ones. 

5.2.1 Momentum integration with insertion operators 

Before we start with specific calculations let us examine how the new factors in the nu- 
merator, originating from the modified Feynman rules of figure 1, affect a general triangle 
integration. Consider a prototype integration of the form: 

d D k (A-k) s 1 1 



(2vr) D (k 2 - ml)" 1 ( (Jfc _ ^ _ ^ ^ _ ag) 2 _ ^ 



V3 



(5.14) 



The initial steps, that is, introduction of Feynman parameters and completion of the square 
in the denominator, may be applied as before and we obtain: 

r(^m) f 1 j_ s u!-i u 2 -i u 3 -i f d D k' (A • (k' + x 2 a 2 + x 3 a 3 )) £ 



dx <5i23 x" 1 x\ 



rwrwrwio^"^ 1 ~ 2 ~ 3 J (2n) D (a/ 2 -A)" 123 

(5.15) 

Before we can proceed, we need to expand the factor in the numerator using binomials: 
TV \ s / \ Jl / • 



^n, t U1 - 1 T n-h+"2-i J2+V3-1 f d D k' (A • k'f 11 
dxd 12sXl x 2 z 3 J {27r) D( kf 2_ A y™ ■ 

Because A 2 = 0, terms involving non-zero powers of (A • k') s ~ n inside the integration 
vanish, and so our expression becomes 

r(,i 23 ) £( s W a2r ; 2(A . a3) ; 2 (5 . 17) 



T{yi) V{u 2 ) T{u z ) ^ Vj 2 



dx (5i23 x 1 1 X 



VI— 1 S-J2+V2-1 ^j2+V3~l 



d D k' 



(2n) D {k' 2 -A) U123 



At this point the regular procedure may be resumed. The end result is simply the normal 
expression for a triangle loop summed over powers of A -a, as in 

where I* c and I cc represent the triangle integration with and without insertion operators, 
respectively. 
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5.2.2 Diagram B-II 

The integral for diagram B-II (figure 3) is given by: 

d D k 2 (A • k 2 f- 1 



(2tt) d [k 2 2 -m 2 Y 4 ((fes+p) 



2 2 



V5 



(5.19) 



(2tt) d ( fc 2^i + k2 f _ m2 y 2 ^ h _ p) : 

The fei integration does not involve operator insertions and so it is performed as usual 
leading to: 



j(_l)^23 r(u 123 -§) [ l dx5 12Z x'i 



a '2__ 2 i3 



(4^)f r(i/ x ) i> 2 ) i> 3 ) 7 (_ X2 (! _ X2 ))^3-f 

/• d D k 2 (A • k 2 ) N - 1 



(5.20) 



£2 y S2 

Following the steps discussed in the previous subsection, we perform the k 2 integration 
involving the insertion operator and obtain: 

( _l)^34 5+ l ( _ A . p) iV-l F(u 12U5 -D) J^^' 1 ) (5-21) 



(4tt) d (m 2 )" 12345 " 15 I>0 T{u 2 ) I> 3 ) I> 4 ) I> 5 ) \ j 



((1 - x B ) (1 - x 2 ) + x 6 y^ 5 - D 
Applying a change of integration variables and introducing a regulator a we get: 

(-ir i2345+1 (-A • pf' 1 I> 123 45 ~D) ^ (N - 1 



(4vr) D (m 2 )" 12345 - 13 r(i/ 1 )r( I / 2 )r(^)r^4)r( I / 5 ) ^ V i 

x I duu^ 13 ^' 1 (1 - ^r 345 ^ 1 t# _1 (1 - u^ 3 " 1 (5.22) 



X d 



(1 - a U! u 4 ) Ul2 ' M5 - D 

After expanding the denominator and performing the parameter integrations we obtain: 

(_l)"1234S+l(_ A , p )*-l r> 1345 -f) 

(4n) D (m 2 r 2345 " D i> 2 ) i> 3 ) i> 5 ) 

j )T(j + v 13 )T(N-j + vt 5 -l) V ' 

x ai T(i + 1^12345 " £>) T(i " 1/13 + f ) r(i + TV - j + Z/45 - 1) r(j + ^123 - f ) 



i=0 



r(i + i) r(i + i/ 45 ) r(i + 7v + z/ 1234 5 - f - 1) 
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Using eqs. (3.5) and (3.6) it is possible to systematically express the above equation in 
terms of Z-Sums and ultimately polylogarithms for generic values of N, f£, and D. For 
simplicity, we will calculate it for i/g = 1 (I = 1, . . . , 5) and D = 4 — 2e. Apart from the 
prefactor (4 ^ 2 >. we get: 

tvo i a V" 1 ( N ~ l \ T{ - j ± l l T U + 1 + £) r(Ar - j) 

( } U\ J J r(j + 2)T(N-j + l) 



oo 



^ r(z + i)r(i + 2)r(i + iv + 2 + e ) ■ lb ' 24j 



At this point the Z-Sum reduction algorithms may be applied, similarly to what was done 
for the cases without operator insertions. We obtain an expression in terms of polyloga- 
rithms, apply the limit a — >■ 1 to remove the regulator, and obtain the final results, which, 
for this particular case, no longer depend on polylogarithms. This calculation has been 
performed for N equal to 2 through 5 in order to reproduce previous results obtained by 
other methods [8], and are presented in table 8. 

5.2.3 Diagram A- VI 

Diagram A- VI involves the insertion of a vertex operator, and so it is a bit more complicated 
than the previous example. We start the calculation with: 



d D k 1 (A • 

(2vr) D ( fc 2 _ m 2)"i ( (fel + k2 fy ^ _ ~f _ m 2y 3 

In this case both integrations include a factor in the numerator. After the k\ integration 
we obtain: 



D 



(4yr) f r^or^r^) 



D\ N-2 n / - n 
— n „■„— n VJ2/ 



jl=0 j2=0 



1 cfa; 5i23 Xj 1 X 



"1 — 1 ^.^2 — 1 ™^3 — 1 



2 



-z 2 (l-x 2 ))" 123 " 



_D 

" 2 



(5.26) 



d D fc 2 (A • k 2 ) N - n+j2 - 2 



. • , ■■> n\Vi I \ 2 o \ J / / , X3p ^ _ m£ 

£2 



(A; 2 - m 2 )" 4 ((fc + p) 2 - m 2 ) 5 ( (k 2 + j 
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We proceed with the fc 2 integration and change parameter variables: 

(_l)^ 23 45+i (A-p) N ~ 2 T(u 12M5 -D) 

(4tt) D {m?) ul2 ^'- D V[yx) I> 2 ) 1> 3 ) 1> 4 ) 1> 5 ) 

N-2 ji , . v N-h+j 2 -2 , . . 

* E E (£) E " a (5 ' 27) 

ii=oj 2 =o V2/ j 3 =o v J6 / 



M ^4-l ( 1 _ U3 )AT-i 1 +i 2 -i 3 + i , 5 -3 ^- J1+ i 2 ^3+^5-3 ^^p+^-f-l 

(l-a«i u 4 )' /12345 " D 
After denominator expansion and parameter integration we get: 

(-1)^+1 {A . p f-2 i ^ * ^ x 

(4^) D (m 2 )^ 5 -^r(z/ 2 )r(i/ 3 )r(^)^ ^ W l ' 1 

X_ "" /J " 2 /iv-j 1+ j 2 -2\ r(j 1 -j 2 +j 3 +^) r(iY-j 1 +j 2 -i 3 +^-2) 



x E 



J3 / r(ji-J 2 +j3 + ^13) r(iV-ji+j2-j3 + l / 45-2) 



i3=o 

x r(i+z. 1234 5-^) r(i+ii-i 2 -^+#) r(i 2 +z/ 24 5-f ) 



E 



x > a 



=0 



r(»+i) r(f+ii+i/4B) 



, r(i+Ar-j 1 +j 2 -j 3 +t/45-2)r(j3+i/ 12 3-f) 

r(i + ^-ji+i 2 + l/i 2345 -2) 

For the purpose of this example we set f£ = 1 (£ = 1, . . . , 5) and D = 4. Apart from a 
prefactor we get: 

V V f ^y 2-2 f « " h + ia - 2\ rOa + ^rQg + ^rQ-x-ja+jg + i) 

h \ h 1 nh-h+h+2) 

, r(jv-ji+j 2 -j3-i) y j £(j + ji - j 2 + 1) r(z + jv - ji + j 2 - j 3 ) , , 
r(Af-ji+j 2 -j 3 ) r(i + i 1 + 2)r(i + 7v-i 1 +j 2 + i) ' l " J 

At this point the Z-Sum reduction algorithms may be applied, similarly to what was done 
previously. We obtain an expression in terms of poly logarithms, apply the limit a — > 1 to 
remove the regulator, and obtain the final results, which, for this particular case, no longer 
depend on polylogarithms. Results for N equal to 2 through 5 are presented in table 8. 

5.2.4 Remaining diagrams 

Having explicitly discussed two calculations as examples, we end this section by presenting 
results for the remaining operator-insertion diagrams required in the calculation of the 
process. Calculations for all diagrams in figure 3 have been reproduced including the first 
four Mellin moments and are presented in table 8. In order to match the previous results 
[8], we did not include a constant prefactor, similarly to what was done to obtain eqs. 
(5.24) (diagram B-II) and (5.29) (diagrams A- VI). 
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N 


2 


3 


4 


5 


A-Ia 


1 

2 


67 
216 


31 
144 


2161 
13500 


A-Ib 


13 

~144 


19 

~432 


17 

~675 


431 
27000 


B-II 


-i + i + T* 
is 4 2 


11 23 117b 


5 11 57 B 
~48e + 96 + ~24~ 


137 949 1377 B 


72e + 144 + 36 


1800e ' 10800 ' 900 


A-III 





_^ + ± + 2£ 

24e 48 12 





90e 270 45 


A-IV 


1 





31 
72 





B-V 


1 1 

"27 + 2 +7B 





5 11 5j E 
~24e + 48 + 12 





A- VI 


1 


1 


65 
72 


29 
36 


B-VII 


1 1 

"27 + 2 +7B 


_l + I + 2£ 

4e 4 2 


5 29 57b 
~24e + 144 + ~Y2 


5 , 7 57J5 
36e 48 18 



Table 8. Results for the first four Mellin momenta for the diagrams shown in figure 3. 



While we will not be presenting results for diagrams involving bubbles as building 
blocks (figure 4), because of their simplicity, we note that the method can be successfully 
applied in all cases. 

6 Conclusions 

We have developed a study of different methods to systematically perform expansions of 
multi-loop calculations leading to concatenated sums, with the intent of applying the Z- 
Sum algorithms in order to obtain final results in terms of multiple polylogarithms. These 
methods involved the use of Taylor expansions or inverse Mellin-Barnes transformations, 
applied either at the level of momentum integration or parameter integration. While most 
of the focus in the literature in recent years has been on the application of approach C 
(see table 1), based on inverse Mellin-Barnes transformations, to specific calculations, we 
decided to cover and develop a way to systematically use Taylor series expansions to obtain 
a general parameter integration. This allowed us to clearly understand the structure of 
the summations obtained, and thus evaluate how often the method can be successfully 
applied, while at the same time making progress towards a generalization of the method 
by identifying missing algorithms. 

We proceeded by surveying the applicability of the Taylor expansion method to one- 
loop diagrams in QFT and found that not all calculations could be performed using the 
Z-Sum machinery as it is, specially if several external legs are off-shell. Since some of 
these diagrams are used as building blocks for multi-loop calculations this demonstrates 
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the limitations of the procedure, with similar difficulties to be expected for box diagrams. 
While this is a concern if the ultimate goal is to develop an automatic computational 
package capable of performing higher- loop calculations, it does not necessarily put this 
approach in a worse position than the alternatives, given that such a general survey has 
not yet been done for the other methods and so their limitations are still unknown. 

Since one would like to explore the possibility of expressing results in terms of multiple 
polylogarithms, we see two different approaches in proceeding with this line of research. 
One would be to try to extend the Z-Sum machinery in order to obtain the necessary steps 
for an entirely systematic reduction using approach A, finding solutions to eqs. (3.19-3.22) 
and their generalizations. The other option would be to forgo the use of Taylor series and, 
similarly to what was done in this work, develop a systematic procedure for performing 
a Mellin-Barnes expansion of a general loop integration using approaches B or C. This 
would allow us to understand the general form of the concatenated sums obtained from 
such methods, finding out how often they can be successfully applied, and which new 
algorithms are necessary for a general solution. 

Finally, having discovered in which situations approach A can be successfully used, in 
section 5 we conclude by applying it to a problem of physical interest, namely the calcu- 
lation of heavy-flavor effects in the structure functions of deep inelastic scattering, which 
involves the calculation of two-loop diagrams with massive internal particles and operator 
insertions. This is an important result since for the first time multi-loop calculations involv- 
ing massive internal particles have been performed using the Z-Sum reduction algorithms 
and Taylor expansions. By doing so, we were able to reproduce results obtained using 
inverse Mellin-Barnes transformations, as presented in [8]. 

This calculation shows that there are in fact cases of multi-loop diagrams involving 
massive particles where the method may be successfully applied. We believe this work 
sets a useful foundation for further investigation into applications of Z-Sum algorithms for 
calculations of loop integrals. 
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A Notation and conventions 



Boolean step function Whenever we use the step function 9 it will have a boolean 
argument as in: 

_ , . I 1 if x = true , . . 

" W = (oifx = Mse ' 

Floor and ceiling function The floor function [^J evaluates as the largest integer 
smaller or equal to x. Similarly, the ceiling function \x~\ equals the smallest integer larger 
or equal to x. Namely: 

[x\ = n, n el, < x - n < 1 , (A. 2) 

\x~\ = n, n € I, < n — x < 1 



B Summation splitting and reordering 

When dealing with concatenated sums, it is useful and sometimes even necessary to shift 
the summation variable and its limits, or reorder sums. An example of such a requirement 
is when one needs to compare two different functions by expressing them in the same basis 
of monomials. In this appendix we list some useful identities that we have repeatedly used 
in obtaining results presented in sections 4 and 5. 

B.l Basic identities with a single sum 

We start with very basic ones involving one sum only: 

n n 

£/(•)=£/(»-<) , (B.l) 

i=0 i=0 
n+i i n 

E = E /(•?') ' (R2) 

j=0 j=0 j=0 

n—i n i 

E^') = E^) + E^' + n - i )-^ n - i ) > (B.3) 

j=0 j=0 j=0 

where f(i) represents any function of i, and n is an integer. 

In eq. (B.3) we added and subtracted terms originally outside the range of the original 
summation, which might lead to badly defined terms (for example involving T(— 1)). These 
will eventually cancel out but a regularizer should be used so that intermediate steps are 
well defined. 

It is possible to express a sum in terms of other summations with multiplied upper 
limit. A trivial example comes from splitting into odd and even sums: 

LfJ LVJ 
E f(i) = E /( 2 *) + E f( 2i + !) > (B-4) 

i=0 i=0 i=0 
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with the inverted identity given by: 



n 1 / 2n 2n \ 

E /(O = o E ^/ 2 ) + E /(</ 2 ) . (B-5) 

i=l \i=l i=l / 



where odd terms simply cancel out. 

These equations can be generalized to: 



1 1 5=2 1 

E/«=E E w+p) > (B-6) 

i=0 p=0 i=0 

and 

E/« = -EEw/(^), (b.7) 

i=l ^ p=0 i=l 

where p and ^ are integers and the coefficient rq is called root of unity [27], and is defined 
by: 

r; = exp(^) , (B.8) 

with the properties: 

ft) j+q = (r p q Y , (B.9) 

and 

q if m mod g = , 
if m mod q ^ . 



9-1 



_ ? limmoa ? = u > (Bin) 

p=0 



B.2 Splitting identities 

The multiplication of Z-Sums is a result of sum splitting: 

n n n i— 1 n J— 1 n 

E E 3) = E E /(*• j) + E E /(*. j) + E /(*> <) • 

i=l 3=1 i=l 3=1 3=1 * = 1 *=1 

A minor limit modification (useful for sums originating from denominator Taylor series) 
gives: 

n n n i n j n 

E E /& j)=EE + E E j) - E /(*. o • ( B - 12 ) 

i=0 j=0 i=0 3=0 3=0 i=0 i=0 

Concatenated sums involving floor function in the upper limit may be split into odd 
and even parts according to: 

a L^J+c LfJ »+Lfj+c L^J i+L^J+c 

E E /&J) = E E /(M.i)+E E + » (B.13) 

t=0 3=0 i=Q 3=0 i=0 3=0 

where a, 6, and c are integer numbers. 
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B.3 Shifting identities 

Other identities obtained from sum reordering are: 



i=0 j=0 


n n—t n n—t 

= EE-^+^) = EE^+^) . 

i=0 j=0 i=0 j=0 


(B.14) 


n n—i 

EE /(*>■?') 


n n—i n i 
i=0 j=0 i=0 j=0 


(B.15) 


oo i 
i=0 j=0 


OO OO oo oo 

= EE^+w) = EE^+^) > 

i=0 j=0 i=0 j=0 


(B.16) 


OO oo 

EE/(^') 

i=0 j=0 


oo oo oo i 

= EE/(j;i) = EE/(i-j.j") • 

i=0 j=0 i=0 j=0 


(B.17) 



Equation (B.15) is simply eq. (B.14) rewritten with f(x,y) — > f(x — y,y), while the last 
two are obtained from the first ones after taking the limit n — > oo. 
Similarly to eqs. (B.16) and (B.17), further shifting leads to: 

oo t oo 

EE/M=EE/(i-j^)=EE/(i-(^-i)i,j) . cb.18) 

i=0 j=0 i=0 j=Q i=0 j=0 

and 

ooi oo 2i / .\ oomi^ / . 1 > \ 

EE/^i)=EE/0+L|j ) i)=EE/0-(i- 1 )-L^J ) iJ • (B.i9) 

i=0 j=0 i=0 i=0 V 7 i=0 j=0 V ' 

All infinite summations being reordered must be absolutely convergent for these identities 
to be valid. 

B.4 Reordering of concatenated sums 

Now we present the algorithm for reordering two sums where the inner sum upper and/or 
lower limits depend linearly on the outer sum's variable. Since this identity does not change 
the arguments within the sum we omit f(i,j)- Expressions involving more than two sums 
can be reordered iteratively, two sums at a time. 
The initial sum is given by: 

b ei+f 

E E > ( B - 2 °) 

i=a j=ci+d 

where a, b, c, d, e, and / are real numbers. We can express the limits of this summation 
as: 

a < i < b , (B.21) 
ci + d<j<ei + f , 
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where i and j are integers. We would like to invert the condition (B.21) and express the 
summation limits in the form: 



fi(a,b,c,d,e,f)<j<f 2 (a,b,c,d,e,f) , (B.22) 
h{j, a, b, c, d,e,f) < i < f 4 {j, a, b, c, d, e, /) , 

which will be rewritten as 

f 2 (a,b,c,d,e,f) f4(j,a,b,c,d,e,f) 

E E • (B.23) 

j=fi{a,b,c,d,ej) i=f 3 (j,a,b,c,d,ej) 

In case we obtain a block of the form (B.22) involving < or > instead of < or > we 
make the substitutions a < i — > [a + lj < i and i < a — > i < [a — 1], which is correct 
whether a is an integer or not. In all other cases when the summation limits are not 
integers the ceiling function for the lower limit and the floor function for the upper limit 
are implied. 

We separate the solution in eight cases depending on the values of c and e. All cases 
will have constant blocks for all values of a, b, c, d, e, and / and also conditional blocks 
with a boolean step function. 

Case 1 e = 0Ac>0 

b f f b cb+d *=r f *¥■ 

E E = E E+ ( c6 +^/) E Y.+^ cb+d> f) E E > ( B - 24 ) 

i=a j=ci+d j=\cb+d+l\ i=a j=ca+di=a j=ca+di=a 

Case 2 e = 0Ac<0 

b f f b \ca+d-l\ b f b 

E E = E E+^ ca +^/) E E+^ ca+d> /) E E > ( R25 ) 

i=a j=ci-\-d j=ca-\-d i=a j=cb-\-d j— 3-d j=ch4-d j— j—d 

c c 

Case 3 c = 0Ae>0 

b ei+f \ea+f-i] b eb+f b eb+f b 

EE= E 5>*(ea + /<d)£ E +0(ea + f>d) £ ]T , (B.26) 

i=a j=d j=d i=a j—d A— j—f j=ea+/ j= 3 ~f 

e e 

Case 4 c = 0Ae<0 

6 ei+f [eb+f-1] b ea+f^r ea+f 

EE= E ^2+o(eb+f<d)^2^2+e(eb+f>d) ^2 £; , (B.27) 

i=a j=d j = d i=a j = d i=a j=eb-\-f i=a 
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Case 5 c>OAe>0 



j — d 

b ei+f cb+d ~ET~ ea+f b 

E E = E E+ E E 

i=a j=ci+d j=[ ea +/+ij i= 3^£ j=\cb+d+l\ i=a 

e 

(cb+d ^ eb+f b ^ 

E E+ E E 
j=ca+di=a j= [ea+f+1] i=t-L j 

(ea+f 1 -^ L eb+f b 

E E+ E E 
j=ca+di=a j=[cb+d+l] i=LzL 



(B.28) 



Case 6 c>OAe<0 



ei+f ^Hr cb+d ^T^T b [eb+f-1] b 

E E - E E+ E E+ E E+ E E 

i=a j=ci+d j=eb+f i=a j = ^ cf-ed ^ %= a j= [cb+d+1] i=a j = ^ cf-ed ^ jj i= a 

( [eb+f-1] *=± ea+f 

E E+ E E 

j=ca+d i=a j=[cb+d+l] i=a j 

i , , 3-d , e 3-f 

cb+d c ea+f e 



3-f ef-ed 
e c — e 



+ 6(cb + d>eb + f) 



(B.29) 



+ 9(cb + d<eb + f)\ E E+ E E 

j=ca+d i=a j=eb+f i=a 



Case 7 c<OAe>0 



ei+f ea+f 



cf—ed 



!l b \ca+d-l] b 



'-f — ed 



11 



EE=EE+EE+E E+ E E 

i=a j=ci+d j = cf-ed i=g j=ca+d i=a j = cf-ed j— j-f j= \ ea +f+] j j = 3~ d 

J c — e ** c— e e r 

(eb+f b [ca+d-1] b 

E E + E E 
j=[ea+f+H i= Lzl j=cb+d i= 3^d 



(B.30) 



+ 9(ca + d> ea + f) 



eb+f b ea+f b 

E E + E E 

j=c.a+d;— 3-f j=cb+d j— LzA 



\ 
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Case 8 c<OAe<0 



b ei+f 



\ca+d-l\ 3 -f \eb+f-l~] b 



E E 



= E E+ E E 



i=a j=ci+d 



j=eb+f j = i~ d j=ca+d i=a 




j=ca+di=a j=cb+d j = 3~ d 




(B.31) 




j=eb+f i=a j=cb+d j— i- d 
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